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ABSTRACT 

We present the first results of our analyses of selected RRLyrae stars for which data have been 
obtained by the Kepler Mission. As expected, we find a significant fraction of the RRab stars to show 
the Blazhko effect, a still unexplained phenomenon that manifests itself as periodic amplitude and 
phase modulations of the light curve, on time scales of typically tens to hundreds of days. The long 
time span of the Kepler Mission of 3.5 yrs, and the unprecedentedly high precision of its data provide 
a unique opportunity for the study of RRLyrae stars. Using data of a modulated star observed in 
the first roll as a showcase, we discuss the data, our analyses, findings, and their implications for our 
understanding of RRLyrae stars and the Blazhko effect. With at least 40% of the RR Lyrae stars 
in our sample showing modulation, we confirm the high incidence rate that was only found in recent 
high-precision studies. Moreover, we report the occurrence of additional frequencies, beyond the main 
pulsation mode and its modulation components. Their half-integer ratio to the main frequency is 
reminiscent of a period doubling effect caused by resonances, observed for the first time in RRLyrae 
stars. 

Subject headings: stars: oscillations — stars: variables: other — stars: individual: KIC 5559631 
(V783Cyg), KIC 3733346 (NRLyr) and KIC 7198959 (RRLyr) 



1. INTRODUCTION 

RRLyrae stars are low-mass stars that have evolved 
away from the main sequence and are burning Helium 
in their core. Their evolutionary stage makes them use- 
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ful tracers of galactic evolution. Like the Cepheids, they 
obey a period-luminosity-color relation and are used as 
distance indicators. RRLyrae stars have typical peri- 
ods of ^0.2 to ~1 d, amplitudes in the optical of 0.3 
up to 2 mag, and spectral types of A2 to F6. Most 
RRLyrae stars pulsate in the radial fundamental mode 
(RRab stars), the radial first overtone (RRc stars) and, 
in some cases, in both modes simultaneously (RRd stars). 
A few RR Lyrae stars are suspecte d to be pulsating in 
highe r-order radial overtone modes (Olech fc Moskalik | 



20091 



1.1. The Blazhko effect 

A large fraction of the RRLyrae stars show a nearly 
periodic modulation of their light curve amplitudes and 
phases on time scales of typically tens to hundr eds of 



days. This so-called Blazhko effect ( Blazhko"|T907 ) is one 
of the most stubborn unsolved problems of the theory of 
radial stellar pulsations. Since it was discovered more 
than a century ago, many hypotheses have been pro- 
posed to explain the modulation. Over the past decade, 
focus has predominantly been on resonance models and 



magnetic models . The oft-quoted resonance models (Van 
'Hoolst, Dziembowski & Kawaler 1998; Dzicmbowski &: 



jMizerski ,2004^ involve nonradial modes of low spherical 
degree (niost likely £ = 1), whereas the magnetic model 



( Shibahashi 
netic held (o: 



2000) proposes that a strong dipole mag- 
the order of 1 kG) inclined to the rotation 



axis deforms the main radial mode to have an additional 
nonradial quadrupole (^=2) component aligned with the 
magnetic field, thus generating the amplitude and phase 
modulation with the rotat i on pe riod 



Rec e ntly, Chadid et al. ( 2004 ) and Kolenberg & Bag- 



nulo (2009) disproved the presence of a strong mag- 



netic field with dipole- like geometry in RRLyr, the pro- 
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totype of the class, and for a sample of RRLyrae stars 
respectively. This implies that the magnetic models in 
their pres e nt for m need adaptation, or are not correct. 



Stothers (2006) proposed a scenario that does not in- 
voive nonradial modes, but convective/magnetic inter- 
action as a cause for modulation. The geometry of the 
magnetic f ield might be too comple x to be detected. Im 



portantly, Jurcsik et al. (2009a I showed that modu- 
lated stars change their mean global physical parameters 
(mean radius, luminosity and surface effective temper- 
ature) over the Blazhko cycle, and phase modulation of 
the pulsation is interpreted to be a consequence of period 
changes. 

Thus far, none of the proposed models successfully 
accounts for all of the features observed in modulated 
RRLyrae stars. The inspiration for new types of mod- 
els should come from the observational side. The Kepler 
mission will be a prime driver for this. 

1.2. Targets 

To select the RR Lyrae targets for the Kepler Mission, 
we searched for all known RR Lyrae stars or possible can- 
didates in the Kepler field. The list was compiled using 
all available variability information contai ned in the up- 
dated catalogs: GCV S ( jSamus et al. 120021), AS AS N orth 
(Pigul skietal. |2009D , ROT SE (lAkerlof et al. |2000D, and 
IMT ( iHartman et af |2004[ ). KlC-lO (KIC = Kepler In- 
put Catalog) was used to get Kepler IDs, coordinates and 
magnitudes. The period range 0.15 — 1.0 d was searched 
in these catalogs . Light curve shapes and the log P vs. 
J — H diagram (Pojmahski & Maciejewski 2004) were 
also utilized for further selection. Finally^ stars with 
close (bright) companions were excluded. RRLyr itself, 
the prototype of the class, is also located in the Kepler 
field. Despite coordinated ground-based efforts we were 
not able to identify any Blazhko star, besides RRLyr it- 
self, in the Kepler field prior to the launch of the mission. 

The Kepler magnitude Kp (wide passband between 
4300 — 9000 A) of the proposed targets lies in the range 
7.9 < Kp < 17.4. A total of 57 RRLyrae candidates, 
including some with uncertain classification, were pro- 
posed for observation. Of these, 48 stars were proposed 
in long cadence (one measurement every 30 minutes) for 
90 d and nine stars in short cadence (one measurement 
every minute) for 30 d. A large fraction (23 out of 51) 
of the stars for which data were released at the time of 
writing turned out not to be RRLyrae stars, but, e.g., 
eclipsing binaries and ellipsoidal variables. 

About 60% of the survey targets will be proposed for 
continued observations after the survey period. 

2. OBSERVATIONS 

The asteroseismic data released to the KASC {Kepler 
Asteroseismic Science Consortium) at the time of this 
writing are long cadence observations. For stars with 
pulsation cycles of typically half a day this sampling is 
sufficient to obtain a good coverage of the light curve, 
but rapid changes and "glitches" in the light curve - as 
have been reported to occur in some RRLyrae stars - 
are missed by this sampling. All of the RRLyrae stars 
with Kepler data have so far turned out to be pulsating 
predominantly in the radial fundamental mode (RRab). 

Note that at the time of writing the pipeline for reduc- 
ing the Kepler data is still being fine-tuned and tested. 



As a consequence, some low-frequency, low amplitude 
variability may be attributable to instrumental effects 
that will be overcome further into the mission. The re- 
sults presented in this letter are valid within this minor 
constraint. 

3. AN EXAMPLE: A BLAZHKO STAR OBSERVED WITH 
KEPLER 

To illustrate the potential of Kepler photometry for 
RRLyrae stars, we present the first frequency analy- 
sis resuhs for the Blazhko star KIC 5559631 (V783Cyg, 
see Fig. 1). Before the Kepler observations the star 



was known to be an RRab star ( Samus et al. 2002 ) 
with a mean Kepler magnitude ot tok = 14.643 and a 
period Pq — 0.620699 d, corresponding to a frequency 
/o = 1.611087 d^^ Blazhko modulation of the fight 
curve was previously unknown. This star was observed 
in long cadence during the first roll of the Kepler sur- 
vey phase between HJD 2454964.0109 (2009 May 12) 
and HJD 2454997.4812 (2009 June 14). (The spacecraft 
"rolls" every three months to allow for continuous illumi- 
nation of Kepler's solar arrays. The first such roll, ending 
the first quarter, occurred on June 18, 2009.) Over the 
33.5-d continuous run 1628 data points were obtained 
for the star. As a consequence of the sampling rate, the 
Nyquist frequency lies at 24.5 d~^. 

A small linear trend was removed from the raw data 
and they were converted to the magnitude scale. We used 
several fr equency analysis soft ware pack ages, inclu ding 
Period04 (Lenz & Breger |2005|), S igSpec ([Reegen | 2007i) 
CLEAN (|Roberts et al. 
„1987t . 



TiFrAn (Csubry 2002) 

and PDM (SteUingwerf , 

After prewhitenmg with the main frequency /o = 
1.61117 ± 0.00003 d~^ and its significant harmonics (up 
to 25/o; although this is past the Nyquist frequency, it is 
necessary, as the higher harmonics are detected by their 
"reflection" into the range below the Nyquist frequency 
by beating with the sampling frequency), we clearly 
see the triplet frequencies that are a typical feature of 
Blazhko RRLyrae stars (see Fig. 1, right panel). The 
first side peak appears at fo + Jb = 1-64735 d~^. This 
yields an initial value for the Blazhko frequency of Jb = 
0.0362 ±0.0006 d"^ or a Blazhko period Pb = 27.6 ±0.5 
d. Subsequent prewhitening reveals additional triplet fre- 
quencies equidistantly spaced around the main frequency 
and its harmonics. The triplet components are detected 
successively and significant up to the 23'''^ order. The 
Blazhko frequency itself is also detected directly from 
our data of KIC 5559631, though the peak in the Fourier 
spectrum at its location is very wide due to the short 
time base of the data. Its detection implies that there is 
a variation of the mean brightness of the star with the 
Blazhko cycle. Though modulated RRLyrae stars are 
commonly fitted with equidistant multiplet structures, it 
is worthwhile to test whether there are departures from 
equidistance. This will be done in our future analyses 
with a longer time base and hence better frequency res- 
olution. 

We fitted the fight curve of KIC 5559631 with the 
main frequency, its harmonics, the significant equidis- 
tant triplets, and the Blazhko frequency. The presence 
of further significant frequencies will be determined on 
the basis of future data of this target. 

Fig. 2 shows the light curve modulation over the 
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Fig. 1. — Left: Light curve of KIC 5559631 folded with the main period of 0.620699 d. The spread in the width of the Hght curve is not 
noise (which is much smaller than the size of the data points), but shows the amplitude and phase modulation over the 27.6-d Blazhko 
cycle. Middle: Fourier Transform of the raw data converted to the magnitude scale. The insert shows the noise level. Right: Fourier 
Transform after prowhitening with the main frequency and its harmonics. The insert is a zoom, on the same scale, around the position of 
the main frequency and its harmonics. Triplet components are clearly detected. 




- 0.6210 



- 0.6205 



0.6200 



T3 
O 



(V 



54965 54970 54975 54980 54985 

HJD-2400000 



54990 



54995 



Fig. 2. — Relationship between the period change and amplitude modulation over the Blazhko cycle for KIC 5559631. The red dots 
show the Kepler data, the red line our multi-frequency fit to the data. The blue line shows an instantaneous period determined using the 
analytical signal method (units to the right side of the panel). The period change cannot be a Doppler shift caused by a companion, as 
that would imply an orbital velocity of the RR Lyrae star about the barycentre of about 200kms~^, requiring an impossibly large mass 
for the companion. 



Blazhko cycle for KIC 5559631. An instantaneous pe- 
riod was det ermined continuousl y by the analytical sig- 
nal method ( KoUath et al. [2002 ) , shown by the blue line 



in Fig. 2. The period change is non-sinusoidal and there 
is a considerable phase lag between the period and ampli- 
tude modulation. Similar behaviour was found by Szabo 



et al. (2009) in the analysis of the CoRoT RR Lyrae 
Blazhko stars. This frequency variability will be fully 
characterized by the long-term Kepler observations and 
will provide an additional new constraint on seismic mod- 
els of RR Lyrae stars, since it must be the result of global 
changes in the star over the Blazhko cycle. 

4. RESULTS FOR THE OTHER RR LYRAE STARS 

For the 28 RRab stars that have been observed by Ke- 
pler so far, we find periods in the range 0.43—0.68 d (main 
frequencies /o between 1.46 and 2.29 d~^), with observed 
amplitudes of the first Fourier component between 0.18 
and 0.44 mag. 

Fig. 3 shows example light curves of non-modulated 
(KIC 3733346, NRLyr) and modulated (KIC 7198959, 
RRLyr itself) RR Lyrae stars in our sample and the 
matching Fourier transforms before and after subtrac- 



tion of the main frequency and its significant harmonics. 
Note the occurrence of modulation peaks and additional 
frequencies in the Fourier spectrum of KIC 7198959. 

Some modulated stars in our sample already clearly 
show higher-order modulation side lobes in their fre- 
quency spectra. Their visibility can be attributed to the 
degree of amplitude/phase modulation (Szeidl & Jurcsik 
2009, Benko et al. 2009) that we will quantify in our 
studies based on larger data sets. 

Besides the modulation components that occur in mul- 
tiplct structures around the main frequency and its har- 
monics, we also see additional frequencies that are not 
at the expected positions of radial modes in RR Lyrae 
stars (see also Section 5.2). They have significant am- 
plitudes surpassing the mmag level (e.g., 8 mmag in 
KIC 7198959). Their nature will be investigated on the 
basis of forthcoming data. 

5. DISCUSSION 

5.1. Modulation statistics 

An important statistic that should help constrain mod- 
els for the Blazhko effect is its occurrence rate. In dif- 
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Fig. 3.— Left: Light curves of two if epkr RR Lyrae stars (KIC 3733346 and KIC 7198959) folded with the main period. Middle: Fourier 
Transform of the raw data converted to the magnitude scale. The insert shows the noise level. Right: Fourier Transform after prewhitening 
with the main frequency and its harmonics. The insert, on the same scale, is a zoom around the position of the main frequency and its 
harmonics. For the modulated star (KIC 7198959) multiplet components are clearly detected. Arrows point at the highest peak connected 
to an additional frequency. 



ferent stellar popul ations the rela tive number of Blazhko 
stars may vary (see Kovacs 2009 and references therein) . 



Previous estimates were that at least 20-30% of the galac 
tic RRab stars and 5-40% of the RRc stars (the high 
RRc incidence rate was fou nd in Omega Cen) are mod 
ulated (Szcidl 1988 

Moskalik & Olech 



2003 



2008 r 



Mizerski 



Moskalik_& Poretti [2003 

Recent results from the 
indicate 



Konkoiy Biazhko Survey (Jurcsik et al. |[2009b[ 
a ~47% occurrence rate ot the Blazhko ettecl (14 out 
of 30 RRab stars). From the excellent CoRoT data an 
even larger fraction of the observed RRLyrae st ars are 
modulated, t h ough the sample size is very small (Szabo 



et al. 2009). Jurcsik et al. (2009b) mention that the 



signihcant increase ot the incidence rate is a consequence 
of the discovery of sm all-amplitude modulation (see, e.g. 



Jurcsik et al. 2006 1 that would not have been detected 



m any previous survey. 

In our sample of Kepler RRLyrae targets, 11 out of 28 
RRab stars are clearly modulated (with large amplitude) 
and for a few more we suspect small amplitude modula- 
tion. Though at this point our numbers are merely in- 
dicative, we obtain an incidence rate of at least 40%, in 
agreeme nt with the findings fr om recent high-precision 



studies (Jurcsik et 
ously thought 



al. 



2009b ), and higher than previ- 

Is every RRab star modulated? What 
fraction of RR Lyrae stars are modulated to within cer- 
tain detection limits? With unprecedented precision and 
the longer Kepler datasets we will investigate these ques- 
tions, too. 

5.2. Additional frequencies 

For a few of our targets, we detect additional frequen- 
cies in our Kepler data that are not at the expected po- 
sitions of the radial modes in RRLyrae stars. However, 
they appear at frequency ratios close to 3/2, 1/2, 5/2, 
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Fig. 4. — Residual spectrum of KIC 7198959 after subtraction of 
the main frequency /o, its harmonics, and the triplet components. 
The dashed lines mark the positions of the half-integer frequency 
components. The highest peak occurs at ~3/2/o. The frequency 
peaks in between the dashed lines are mostly higher-order multi- 
plet components present in the spectrum around /o and its har- 
monics. We did not subtract these yet, as the IBlazhko period 
of KIC 7198959 is only partially covered by these data. The insert 
shows the alternating heights of the maxima around minimum light 
for the star, connected to the additional frequencies. 

7/2, with the radial fundamental mode /o (Fig. 4). These 
ratios are reminiscent of resonance phenomena, specifi.- 
cally period doubling bif urcation as described by Moska] 



lik & Buchler 



( 1990 1 and 



Buchler & Moskalik (1992) for 



Cepheid an d 3LHer models (see also I' lg. 3 m |Moskalik| 



fc Buchler j [l991 
period doub 



for comparison with our Fig. 4). Such 
ling would also result in alternating heights 



of maxima, as we observe in these targets (Fig. 4) . 
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In his recent review of the Blazhko effe ct, and our (lack 



of) understanding of the phenomenon, Kovacs (120091 
stresses that more attention should again be pa id to mod- 
els based on radial mode resonances ( see, e.g. , Goranskij 



etal. 2009). After the models by |Borkowskin ( |1980| ) anJ 
Moskahk (1986[ ) were dismissed by lack of conhrmation 
from hydrodynamical models, most attention has been 
on resonances involving nonradial modes. 

In several recent studies of RR Lyrae stars additional 
frequencies w ere detected, e.g., in phot ometry of the RRd 
star AQLeo ( |Gruberbauer et al. [20"07 ) from the MOST 
satellite, and m the V1127A ql data from the CoRoT 



mission (Chadid et al. 2009). These frequencies were 
not found at positions where radial modes or radial reso- 
nances would occur. They may, therefore, be connected 
to nonradial modes. Nonradial modes play a crucial role 
in se veral m odels for explaining the Blazhko modulation 
(e.g., [Dziembowski fc Mizerski ^^2004) . Their presence 
may, however, be unrelated to tke Blazhko effect. De- 



spite several claims (Olech et al 



[][T999) of detection and attempts |kolenberg ||2002[ )"at 



1999 Chadid et al. 



tricir identification, there has been no unambiguous proof 
of nonradial modes in RR Lyrae stars so far. We hope 
that through the ultra-precise data from the Kepler mis- 
sion we will be able to obtain certainty about their exis- 
tence. Confirmation and exploration of nonradial pulsa- 
tion would be extremely important to shift the current 
radial (1-dimensional) theoretical modeling to a nonra- 
dial (multi-dimensional) hydrodynamical description of 
RR Lyrae pulsations. 



5.3. Quality and advantages of Kepler photometry 

From the first releases of Kepler data it is obvious 
that Kepler photometry is stable enough to be tied to 
ground-based photometry, even if data are not taken si- 
multaneously. Ground-based follow-up photometry and 
spectroscopy of selected Kepler targets is underway to 
improve the information in the Kepler Input Catalogue. 
High-resolution multi-color single-epoch imaging of the 
targets is important to mitigate the 4 arcsec sampling of 
the Kepler photometeT . The latter may cause some prob- 
lems concerning proper identification, determination of 
the total pulsational amplitude, and blending. 

Although the CoRoT satellite is carrying out a similar 
mission to Kepler, and MOST is also delivering valu- 
able results in asteroseismology, Kepler has significant 
advantages for the study of RR Lyrae stars. These in- 
clude the large aperture of the primary mirror: 95-cm 
diameter (vs. 28-cm for CoRoT and 15-cm for MOST) 
allowing the observation of fainter stars and higher preci- 
sion photometry, the potential for obtaining observations 
over a longer time base: 3.5-5 yr (versus a maximum of 
150 d for CoRoT and about 8 weeks for MOST in or- 
der to cover also longer cycles for modulated RR Lyrae 
stars (e.g., RRLyr's 4-yr cycle), and the larger field of 
view: 105 square degrees which allows the observation 
of a significantly large sample of RR Lyrae stars simul- 
taneously. Finally, thanks to Kepler^s Earth-trailing or- 
bit (vs. CoRoT's and MOST's low-Earth orbits), Ke- 
pler is spared from transits through the South Atlantic 
Anomaly (SAA). 



6. SUMMARY 

The first releases of Kepler data have yielded light 
curves with unprecedented quality. In this Letter we 
have attempted to illustrate the tremendous potential 
of Kepler for studies of RR Lyrae stars. 

The major results and prospects for further analysis 
can be summarized as follows: 

• The fraction of modulated stars in our sample of 28 
observed stars so far confirms that the Blazhko ef- 
fect is m uch more co mmon than previously thought 



(see also Jurcsik et al. |2009b| ) . We can already fix 
a lower limit of ~40% tor the incidence rate in our 
sample. The Kepler data will help to determine 
a lower limit for the modulation amplitude and to 
constrain the models. 

• We observe additional frequencies, beyond the 
main frequency and classically expected modula- 
tion components, in the frequency spectra of a 
number of modulated RR Lyrae stars. Their ratio 
to the main frequency suggests resonance effects, 
in particul ar a period doubling bifurca tion as de- 



(1990). It is the 



scribed by Moskalik fc Buchler 
first time such frequencies are reported m RR Lyrae 
stars. 

• The star KIC 5559631 (V783 Cyg) presented in this 
Letter has its modulation cycle (27.6 d) covered by 
the Ql data. For the other Blazhko stars the mod- 
ulation periods are longer. Forthcoming data with 
a longer time base will allow for a multitude of 
analyses: besides the classical ways of analyzing 
Blazhko stars (Fourier techniques, phase disper- 
sion minimization), we will investigate the relative 
contribution of p hase and amplitude modulation in 
RRLyrae stars (Szeidl & Jurcsik 2009 Benko et 
al. |[2009l ). 



• In our forthcoming analyses, we will search for sec- 
ond and higher overtone (radial) modes in the time- 
series photometry, as well as potential nonradial 
modes. 



Several authors (|LaCluyze"||2004| iKolenberg et al 



2006 



Jurcsik et al. 
1999 Goranskij et al 



.yze ZUU4UJ 
feOO^bilC lc: 



ylcment lioran 
12009) have found 



that the modulation over the Blazhko cycle can be 
very unstable. The long time base of the Kepler 
data will allow us to study the stability of pulsa- 
tion and modulation with unprecedented precision. 
The Kepler data also have the required precision 
to investigate irregularities in the pulsation of non- 
modulated RRLyrae stars. 

• We will also investigate features connected to shock 
wave phenomena, such as the change in slope at 
mid-rising light and the am plit udes and durations 
of the hump and bump (see Gillet fc Crowe ,|1988| 
in both stable RRab stars and during Blazhko cy- 
cles. 

• Theoretically predicted strange mode pulsation 
( Buchler fc KoUath [2001 ) , which has not been ob- 
served in RRLyrae stars yet, may also be present 
in the data. 
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• Further possible applications of the 3.5-5 yr Ke- 
pler data are the measurement the secular variation 
of the pulsation period(s), directly li nked with the 



evolut ionary status of a star (see, e.g., Poretti et al 
"] |2008[ ), and the detection of substellar companions 



with I he timing method, as was done with another 



class of hor izontal branch pulsators (e.g., Silvotti 
et al. |[2007l). 
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